Effect of Symmetry Breaking Perturbations in the One-Dimensional SU(4) 

Spin-Orbital Model 
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We study the effect of symmetry breaking perturbations in the one-dimensional SU(4) spin-orbital 
model. We allow the exchange in spin (Ji) and orbital (J2) channel to be different and thus reduce 
the symmetry to SU(2) g) SU(2). A ma™*tic field h along the 5*^ direction is also applied. Using 
the formalism developped by Azaria et altJ we extend their analysis of the isotropic Ji = J2, h — 
case and obtain the low-energy effective theory near the SU(4) point in the generic case Ji 7^ J2, 
h ^ 0. In zero magnetic field, we retrieve the same qualitative low-energy physics than in the 
isotropic case. In particular, the isotropic massless behavior found on the line Ji = J2 < K/A 
extends in a large anisotropic region. We discover however that the anisotropy plays its trick in 
allowing non trivial scaling behaviors of the physical quantities. For example the mass gap M has 
two different scaling behaviors depending on the anisotropy. In addition, we show that in some 
regions, the anisotropy is responsible for anomalous finite size effects and may change qualitatively 
the shape of the computed critical line in a finite system. When a magnetic field is present the 
effect of the anisotropy is striking. In addition to the usual commensurate-incommensurate phase 
transition that occurs in the spin sector of the theory, we find that the field may induce a second 
transition of the KT type in the remaining degrees of freedom to which it does not couple directly. 
In this sector, we find that the effective theory is that of an S0(4) Gross-Neveu model with an 
fe-dependent coupling that may change its sign as h varies. 

PACS No: 75.10.Jm, 75.40.Gb 



I. INTRODUCTION 

In the very past years, there has been an intense in- 
terest devoted to one dimensional spin-orbital models 
[n|. The main reason stems from the recent discovery 
of the new quasi one dimensional spin gapped materi- 
als NaaTiaSbaO |] and NaVaOg [|. It is believed that 
the unusual magnetic properties observed in these com- 
pounds should be explained by a simple two-band Hub- 
bard models at quarter filling. At this filling, and in 
the large Coulomb repulsion the effective Hamiltonian 
simplifies greatly and is equivalent to a modcj_ of two in- 
teracting spin one-half Heisenberg models 



i 

+ k(s,-S,+i]{T,-T,+i] (1) 



where Si and Ti are spin-I/2 operators that represent 
spin and orbital degrees of freedom at each site i and 
Ji(2) and K are positive. It is important to notice that 
(||) is not only SU(2) invariant in S space but also in T 
space. For generic couplings ^1(2) and K, the Hamilto- 
nian (m is SU(2)s(g) SU(2)i symmetric. Two particular 
cases are of interest. First, when Ji = J2 there is an 
additional Z2 symmetry in the exchange between Si and 
Ti. The other important case is when Ji = J2 = K/A 
in which case the Hamiltonian is SU(4) invariant [p|. In 



fact, the Hamiltonian dll) is a simplified version of the 
most general case. Indeed, depending on the microscopic 
couplings, there can be other terms that breaks the SU(2) 
symmetry in both spin and orbital sectors. The choice 
of studying (^ is dictated by the fact that since it re- 
tains some symmetries it is a simple starting point from 
which one can expect to gain some insight before tackling 
with the more general case. In this respect, the model 
(|l|) describes the simplest physically relevant symmetry 
breaking pattern SU(4) -^ 811(2)^(8) SU(2)t. 

The Hamiltonian dl^) can be interpreted in terms of a 
two-leg spin ladder coupled by a four spins interaction. 
Such an interaction can be generated either by phonons 
or, in the doped state, by conventional Coulomb repul- 
sion between the holes M. Some of the properties of 
(|l|) are well established in the weak coupling limit. In 
the limit K <^ Ji(2), the Hamiltonian (nl) describes a 
non-Haldane spin liquid where magnon excitations are 
incoherent m. 

The strong coupling regime, K ^ Ji{2), has just be- 
gin to be investigated P,p||q-[lO|] . From the theoretical 
point of view the situation is ackward. Indeed, as stated 
above, at the special point Ji ~ J2 — K/A, the Hamil- 
tonian (|l|) has an enlarged SU(4) symmetry Q and is 
exactly solvable by the Bethe-ansatz [O. The model is 
critical with three gapless bosonic modes and flows in 
the infrared towards the Wess-Zumino-Novikov-Witten 
(WZNW) model SU(4)i ||l|l. It is conformaly invariant 
with the central charge c — 3. Therefore, there should be 



a qualitative change in the physical behavior of (m) when 
going from small to large K. ^From the theoretical point 
of view this situation is striking since this means that 
one cannot go continuously from weak to strong coupling. 
This is a manifestation of the Zamolodchikov c-theorem 
|l3| that states that starting at K = with two decou- 
pled S — 1/2 chains with the central charge c — 2 (two 
gapless bosonic modes) one cannot flow, in the Renor- 
malization Group (RG) sense, towards the SU(4) point 
which has c — 3 (three gapless bosonic modes). Very 
recently, a new approach to tackle with the strong cou- 
pling regime has been developped by Azaria et al p4[ . 
The idea is to start from the strong coupling fixed point 
SU(4)i and to perturb around it. 

This strategy has been applied to the symmetric line 
J = Ji = J2. It was shown that when J < K/A a small 
deviation from the SU(4) point is irrelevant and thus the 
low energy physics is governed by the SU(4)i fixed point. 
In contrast, when J > K/A the interaction is relevant 
and a gap opens in the spectrum. The system dimerizes 
with alternating spin and orbital singlets. In addition 
it was argued that the SU(4) symmetry was restored at 
long distance and that the low energy effective Hamil- 
tonian was that of the SO (6) Gross-Neveu (GN) model. 
The low-lying coherent excitations was then shown to 
be fermions that transform as an antisymmetric tensor 
of rank two of SU(4). These excitations are coherent 
with wavevector near 7r/2. 

The purpose of this work is to extend the analysis 
to the asymmetric case (Ji 7^ J2) and to inquire how 
anisotropy modifies the low energy physics described 
above. This paper is organized as follows. In Section 
II we present the tools that are necessary to explore the 
vicinity of the SU(4)i fixed point and derive the effective 
low-energy theory associated with (^. A detailed renor- 
malization group analysis is then presented in Section III. 
We obtain the phase diagram in the the plane (Ji, J2) 
and discuss the asymptotic behaviors of the mass gap M . 
A discussion on cross-over effects linked to the anisotropy 
is also presented. In Section IV we investigate the effect 
of a magnetic field on the anisotropic model. Finally in 
the section V we summerize our results and present some 
technical details relative to the computation of the mass 
gap in the Appendix. 



II. THE LOW ENERGY EFFECTIVE FIELD 
THEORY 

A. The SU(4) Heisenberg chain 

Our approach is very similar to the description of the 
S=l/2 Heisenberg spin chain at low energy from the spin 
sector of the repulsive Hubbard model at half filling [n2| . 
To this end, let us introduce the SU(4) Hubbard model 
with C/ > 0: 



n. 
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2 ^-^, 

iahaa 



^ ^i+laa^io-<^ ~^ "'•"^' 



laan.t,^' {l-5ab5^^'). 



(2) 



Here c\^^ denotes electron creation with channel (orbital) 
a = 1, 2 and spin <t =|, [ at the i-th site. The occupation 
number is defined by: n. 



laa iaa ''^^ 



In the limit of 



large positive U, and at quarter filling, it is not difficult 
to show that (0) reduce to (0) with the identification: 



5^,--^^^ 



2 /^"^Jaaf^a/jCja/J 



Jj — 2 7 ,^'iaoJabCiba 



(3) 



where a (resp. r) are the Pauli matrices acting in the 
spin (resp. orbital) space. The low energy physics can 
be described in terms of right movers Raa and left movers 
Laa fermions which correspond to the lattice fermion ciaa 
in the continuum limit: 
c ■ 

_ ~ Raa {x) exp (ikpx) + Laa {x) CXp {—'ikpx) , X — ittQ 



(4) 

where ag is the lattice spacing and the Fermi wavevector 
is defined by: kp = Tr/AaQ. At this point we bosonize 
and introduce four chiral bosonic fields ^aaF!..L using the 
Abelian bosonization of Dirac fermions |15(| : 



Ran = J^ exp (iVAn^aaB^ 



y/Tnao 



Lr, 



\/2Trao 



exp (-iVAn^acrL] 



(5) 



where the bosonic fields satisfy the following commuta- 
tion relation: 



[^aaR,'^ba'L] = T'^ab'^^^' 



(6) 



so that {Raa{x), Lf^y (y)} = 0. The anticommutation 
between fermions with different spin-channel indexes is 
insurred by the presence of Klein factors (here Majorana 
fermions) Haa with the following anticommutation rule: 



{Kaa,Hba'} = 2(5afc(5g 



(7) 



As in the the solution of the two-channel Kondo effect 
by Abelian bosonization 119] it is suitable to introduce 
the physically transparent basis: 



$/= 2(^11 +*U-*2T -*2i) 

'^sf^^i^n-^n- '^2^ + ^21)- (8) 



In this new basis, the total charge degree of freedom is 
described by the bosonic field <i>c while the other "spin- 
orbital" degrees of freedom, are faithfully bosonized by 
the three bosonic fields <I>s,$/,$s/. It is now straight- 
forward to obtain the continuum limit of the Hubbard 
Hamiltonian (0): 



Ti-u = Tic + 'Hsf 



(9) 



where 



Wc = Y ((^-*=)' + ^d.Qcf) + ^ (a,$c)2 , (10) 



and 



^^/= E {\{{dAf + {d.Qaf 



a=sj,sf 



Uao 2 



cos V 47r<I>, cos V 47r$ 



^f 



U 



+ cos VAn^f cos V47r$s/ + cos V47r$s cos v47r$s/ ) (H) 

As in the SU(2) Heisenberg chain, spin and charge de- 
grees of freedom separate. Notice however that at this 
order in U there are no umklapp terms in the charge 
sector since we are at quarter filling and the Akp contri- 
bution to the effective Hamiltonian oscillates. Umklapp 
terms will arise at higher order in perturbation theory 
and will be responsible of a Mott transition at a finite 
value oiU — Uc jl^. Assuming that [/ >> C/c, we focus 
now on the spin-orbital sector. 

The interaction term in Hsf has scaling dimension 2 
and is therefore marginal. This term is nothing but the 
SU(4) current-current interaction. The 15 SU(4)i cur- 
rents Jj^ and Jl can be expressed in terms of the three 
bosonic fields $s,^ and $s/ and the spin-orbital part of 
the Hamiltonian (O) takes the form: 



n 



2ttv^ 



sf 



where 



J2 iJRjR + J£J£) + '29s E ^R-^L 



a=l,15 



E 

a=1.15 



(12) 



gs = -Van 

Vs — vp — C/ao/27r. (13) 

The first term in Eq. (|l2| ) is just the Sugawara form 
of the SU(4)i WZNW model while the second term is 
the marginal current-current interaction. When U > 0, 
gs < 0, the current-current interaction J^J^ is thus ir- 
relevant, as a consequence the spin orbital sector is de- 
scribed by the SU(4)i WZNW model 111. 



B. Majorana representation and the SO(6) 
Gross-Neveu model 

As first emphasized by Shelton et al |Qq| in their study 
of the two-leg spin 1/2 ladders it is very convenient to 



formulate spin ladders problem in terms of real (Majo- 
rana) fermions. This can be done by refermionising the 
three bosonic fields ^s,*!*/ and <i>s/. Let us introduce the 
six Majorana fermions ^", a = 1..6 as follows: 






ie 



ie 



m 



if 



'R{L} 



'R(L) 



R{L) 



R{L) 



: exp ( ±i-\/47r$ 

: exp (±iV^^fR{L) 

■.CXp(±iViTT^sfR{L 






(14) 



where rji are Klein factors. With all these relations at 
hand, one can rewrite the Hamiltonian (O) in terms of 
six Majorana fermions: 



n 



s/ 



J2(^R^-^R-^i^-^L^+9sT.^^^i (15) 



i<j 



where we have introduced the energy density of the dif- 
ferent Ising models: iea = Ka = ^r^l- The Hamiltonian 
(|l|) is nothing but that of the SO (6) GN model with a 
marginally irrelevant interaction when C/ > 0. Thus in 
the far infrared the six Majorana fermions decouple and 
remain massless. This is the S0(6)i fixed point. 

The above result assumes that U is small and the ques- 
tion that naturally arises is whether it can be extended 
to large values of U where (pi) reduces to (El) . For exactly 
the same reasons as for the SU(2) Heisenberg chain the 
answer is positive. We know from the exact solution that 
the SU(4) model is critical with three massless bosonic 
modes or equivalently six massless Majorana fermions. 
We know from Conformal Field Theory that the fixed 
point Hamiltonian can only be the SU(4)i ~ S0(6)i 
WZNW model. The marginal interaction ^^ 



r\ij !%■} lb 



<j "'*"'j 



the only one that respects both the SO (6) symmetry as 
well as translation invariance. Therefore in the vicinity 
of the fixed point the SU(4) Heisenberg model will be 
given by: 



n 



J2{CRd.CR-eLd.^l) + 2G,J2^.K, (16) 



a=l 



Kj 



where the spin velocity Ug and the coupling G3 < are 
unknown and non universal parameter that could be ex- 
tracted from the exact solution. The only thing that 
happens when going from small U to large f/ is a renor- 
malization of the gs and Vs ■ 

The Majorana description used here is extremely use- 
full to understand the symmetry properties of our model. 
Indeed for example, one can define the spin and orbital 
triplets: 



(tR{L) = (C ;C jC 



RiL) 



'R{L) 



(17) 



These quantities transform like a vector under spin 
S0(3)s and orbital S0(3)t rotations. These correspond 
to the SU(2)s and SU(2)f transformations acting on the 
operators S and T respectively. 

To get a complete description of the SO (6)1 fixed point 
one needs the continuum expressions for the effective spin 
and orbital densities in terms of the Majorana fermions: 

S = JsR + JsL + cxp {mx/2ao)J^s + H.c. + (-1)^/^" hs(x) 
f = JtR + JtL+ cxp (mx/2ao)N't + H.c. + {-lf''° nt(x) 

(18) 

where Js.t is the uniform, A; = 0, part and Ms,t and 
fis^t are the 2kF = 7r/2ao and 4:kF = t^/clq contribu- 
tions. Notice that in contrast with the SU(2) Heisenberg 
chain, the spin density has three oscillating components. 
The reason for this comes from Conformal Field The- 
ory. Indeed the different oscillating components of the 
spin density are primary fields of the SU(4)i WZNW 
model. There are three of them with scaling dimensions 
(3/4,1,3/4) [|l9|. They all belong to the representations 
(building blocks) of SU(4) with Young tableau consisting 
of a (a=l,2,3) boxes and one column. In particular the 
staggered components at 4ki? = tt, n^.^^, are components 
of a antisymmetric tensor of rank 2. 

The uniform components express in terms of the Ma- 
jorana fermions as follows: 



J, 



sfl(L) 



Ji 



tR{L) 



^sR{L) A ^sR{L) 
; £.tR(L) A £,tR{L) 



(19) 



Notice that in contrast with the SU(2) Heisenberg chain, 
the uniform part of the spin densities are SU(2)2 currents. 
The expressions for the 4^^ = tt/oo densities are given 
by 



fis = i-B ^sR A £_sL 
fit = iS itji A itL 



(20) 



where i? is a non universal constant. Their scaling di- 
mension at the SO (6)1 fixed point is A^ = 1. Both den- 
sities (O) and ( pO| ) are rather simple when expressed in 
terms of the Majorana fermions. This is not the case 
with the 2kp — 7r/2ao densities JVs,t that are non local 
in the Majorana fermions ^°'. Indeed they involve order 
and disorder operators aa and fia of the six Ising models 
that are associated with the six Majorana fermions. The 
expressions of both Ms,t are lengthly and we shall give 
here only the z component that will be sufficient for our 
purpose: 

TV/ = A {iHlH2Cr3CriO-50-6 + CriCT2/i3Ai4Ai5M6) 

TV/ = A (io-lO-2/^3M40-50-6 + AilM2Cr30'4Ai5Ai6) (21) 



where A is also a non universal constant. Since at the 
free fermion point, the order and disorder operators have 
scaling dimension 1/8, the 2kp densities TVs,t have the 
scaling dimension A^/2 = 3/4. 

This completes the continuum description at the SU(4) 
point. The theory is critical and flows to the SO (6)1 fixed 
point. There is a marginally irrelevant correction which 
magnitude is nonuniversal and is given by the unknown 
coupling G3 < 0. This is in complete agreement with 
the non-Abelian bosonization approach of Affieck p2[ . 
With the continuum expressions of the spin and orbital 
operators at the SU(4) point in terms of the six Majo- 
rana fermions and the associated Ising models we shall 
be able us to investigate the properties of (|l]) close to 
the SU(4) symmetric model. But before moving to this 
point, let us stress that the effective theory depends on 
three unknown parameters A and B and G3 which can 
be in principle extracted from numerical studies. In the 
following, we shall assume that the above description still 
holds for small deviations of the SU(4) point. The only 
modification being a renormalization of the non universal 
constants A, B and G3. 



C. The SU(2)(g)SU(2) symmetry breaking 
perturbation 

We shall now derive the effective low energy theory 
associated with (^ for arbitrary values of Ji and J2 close 
to the SU(4) invariant point given by Ji = J2 = K/4:. 
To this end, let us parametrize the couplings as follows: 



J2 = J + G2 



(22) 



where both Gi and G2 are much smaller than K. The 
Hamiltonian dl^) can then be written as: 

'^ — '^su(4) + Gi 2_^ Si ■ Si+i. + G2 ^^Ti • Ti+i. (23) 

i i 

Using the low energy description of the spin-orbital op- 
erators (|l^), (|2y) and (pi|), one can expand (|2|) around 
the S0(6)i fixed point: 



n 



iXsR ■ dxisR - 6l ■ dx^LJ + 51 («1 + K2 + H^f 
-Z- iXtR ■ dx^tR - 6l • d-MLJ + 92 (k3 + K4 + K^f 



lUs 

~2' 



+ 2G3 (ki + K2 + Kq) (k3 + K4 + K5) 



(24) 



where 



5i = Gi + G3 
52 = G2 + G3 



(25) 



and the two renormalized velocities Us and ut are given 
by: 



Us 
Ut 



2Gi/tt 
2G2/^. 



(26) 



In the above equations and in the remaining of this paper, 
we include the effect of the 4ki? components of the spin 
and orbital densities in a redefinition of the couplings: 
G(i,2) -^ (l + 5^)G(i.2)- The Hamiltonian @ describes 
two marginally coupled S0(3) Gross-Neveu models: one 
in the spin channel described by the three Majorana ^^ 
and one in the orbital channel described by the three 
Majorana ^(. The situation at hand is to be contrasted 
another time with the one encountered in the study of 
spin ladders. In the latter models deviation from crit- 
icality leads, in general, to relevant perturbation and a 
gap always opens independently of the sign of the cou- 
plings. There are however notable exceptions where frus- 
tration plays its trick. This is the case of some three-leg 
frustrated ladders where for some particular values of the 
couplings only marginal interaction remains in the effec- 
tive action. As a result a non trivial critical state, the 
so-called "chirally stabihzed" hquid, shows off Ig^l- In 
the model the situation is even more striking since we 
find only marginal interactions in a finite region of the 
couplings. This is mainly due to the fact that the frusta- 
tion is maximal in the strong coupling region. A direct 
consequence of the marginality of all the interaction is 
that we expect the phase diagram to result from a deli- 
cate balance between the different terms entering in (24). 
This is why it is now worth discussing the effect of the 
three different terms in (p4|). 

Consider first the case where G3 = 0. Then we are left 
with the two decoupled S0(3)s and S0(3)t GN models 
which are exactly solvable. At issue is the sign of Gi 
and G2. When Gi > and G2 > a gap opens in 
both spin and orbital channels. The spectrum consists of 
kinks and antikinks (there are no fermions) ^l|. When 
Gi < and G2 < 0, the interaction is irrelevant and the 
model flows towards the isotropic S0(6)i ~ SU(4)i fixed 
point with the central charge c — 'i. When G1G2 < 
one of the SO (3) GN model will become massive while 
the other one will flow towards the S0(3)i fixed point. 
The Hamiltonian will remain critical with three massless 
Majorana fermions leaving the whole system with the 
central charge c = 3/2. 

Now the physically relevant question is whether or not 
this scheme survives a small negative G3. Indeed, it 



would be not correct to neglect the last term in (24). 
First of all, from the point of view of the lattice Hamilto- 
nian ( p3| ) , the Hamiltonian (G4) has to be thought as the 
effective Hamiltonian obtained by integrating out high 
energy modes up to a scale where one sits close enough 
to the SO (6)1 fixed point where the continuum limit can 
be taken. Thus, generically G3 is not zero. The second 
reason is that all other interactions are marginal. In such 
a case, it is well known that operators that are naively ir- 
relevant may become dangerous and strongly modify the 
physics in the infrared. Therefore, eventhough G3 < 0, 
one has to keep it and analyze with care (e3) with all cou- 



plings different from zero. As we shall see, the strong ten- 
dency to the SO (3)1 criticality in the regions G1G2 < 
will be spoiled in most of the parameter space. There will 
be though still a finite region where the Hamiltonian will 
be critical but with an approximate SO (6) symmetry. 



III. RENORMALIZATION GROUP ANALYSIS 

A. The renormalization group flo^v 

The RG equations for the couplings entering in (|24 
are given at leading order by: 

91 = 9I + Gl 

TTUs TTUt 

92 = 52 + Gl 

TTUt TTUs 



n Us Ut 



(27) 



In the above equation, G means dG/dt where t = log(A) 
is the RG scale. It is more convenient to express the set 
of Eqs. (p^) in terms of the couplings that enter in the 
lattice Hamiltonian (p3|): 



Gl — Gl — 2G2G3 
G2 = G2 — 2G1G3 
G3 = 2G3(Gi+G2 + 2G3), 



(28) 



where, for simplicity we have made the following redefi- 
nition: 



Gl ^ l/iTy/UsUt {aGi - (1 - q;)G3) 
G2 -^ l/TT^^i:^t{l/aG2 - (1 - l/a)G3) 
G3 -^ l/ny^UsUtGs 



(29) 



with a — y^ut/ug. At the leading order, the velocities Ug 
and Ut also renormalize: 



Ut. 



lis = ~6usG^ I{ — ) 



Ut = -GutGi X{ — ) 



Ut 



where 



I{x) 



1-x 

1 + x' 



(30) 



(31) 



The RG equations (pSl) and (p^ can be exactly solved 
and reduce to a single differential equation. To show this 
let us introduce the following variables: 



/i — G1G2 + G1G3 -|- G2G3 
d = Gl - G2 

s = Gl + G2 + 2G3 



(32) 



Eqs. 



and (30) greatly simplify to 



and 



11 = 113 

d = ds 
G3 = 2G3S 



lis -ut = -GGl {us 
UgUt = constant. 



Ut) 



(33) 



(34) 



While the meaning of both s and d is clear, the physical 
interpretation of the variable /i is not straightforward. As 
we shall see in the following, /i measures the departure 
from criticality. 



The solution of both (^ 
tained: 

f,{t) = M(0)X(t) 

d{t) = d{0)X(t) 

G3it) = G3{0)X^it) 



and (p4) are then easily ob- 



{us - ut){t) = {u, - ut){0) exp ^-6 GJiO) f X\T)dT 

(35) 



{UsUt)(t) = {UsUt){0) 

where 



X{t) 



exp / s(r)dT 
'0 



is the solution of the differential equation: 



= P{X). 



(36) 



(37) 



In the latter equation P{X) is a fourth order polynomial 
that depends on the initial conditions of the flow: 



P{X) = 4:Gl{0)X-^ + d'^{0)X^ + 4pi{0)X. 



(38) 



Once the solution X{t) of (37) is known, the behavior of 
the RG flow is completely determined. In particular the 
time evolution of couplings Gi{t) and G2(t) is given by: 

1 



G(i.2)(i) = - [eis)^PiX{t))-2G3{0)X\t) ± d(0)X{t) 

(39) 

where e(s) = sign(s). 

As in the XXZ model we distinguish three different 
phases A, B and C that are separated by the two surfaces 
defined hy fi = 0. 




t 

FIG. 1. RG flows for the variable X(t) in the different 

phases A, B and C. In the phase C, X(t) reaches a minimum 

value X* at a scale t*. When s(0) < 0, the minimum is 

reached for a physical time f > 0. 



- The A Phase. There one has ^(0) > and s(0) > 
and A"(i) ^ 00 as i ^ 00 (see Fig. 1(a)). All couplings 
are relevant and a gap opens in the spectrum. In addi- 
tion, the velocities Us{t) and ut{t) goes to a fixed value 
u'* = Wf . Looking in more details in the large t behavior 
of the couplings, one finds that Gi{t) ~ G2{t) ~ — 2G3(i) 
so that the effective interaction in the far infrared takes 
the suggestive form: ~ {{ki + K2 + kq) — (k^ + K4 + ^5))^ 
which after the transformation ^sR —* ^sR, ^sL —> —^sL 
acquires a manifestely S0(6) invariant form ~ (^i Ka)^- 
One may therefore be tempted to conclude that the 
SO (6) symmetry is restored in the far infrared and that 
the low energy spectrum of the Hamiltonian would be 
approximatively that of the S0(6) GN model. However, 
as pointed out by Azaria, Lecheminant and Tsvelik [E2| , 
one cannot conclude on the basis of perturbation theory 
alone. In fact looking at the flow in the infrared is not 
sufhcient and can be misslcading. The point is that in the 
A phase, since X(i) ^ as t ^ — 00 all couplings goes to 
zero in the ultraviolet and the theory is asymptotically 
free. From the field theoretical point of view this means 
that a well defined renormalized theory with three dif- 
ferent renormalized couplings Gir,G2r and G^r exists 
p3. Therefore, it is most likely that the S0(6) symme- 
try is not restored in the A phase. The situation here is 
very similar to the XXZ model where in the A phase (the 
Ising region), though there is an apparent SU(2) symme- 
try restoration in the far infrared, the exact solution tells 
us that this is not the case [ p3[ . 

- The B Phase. There one has Ai(0) > and s(0) < 
and the flow in the X variable is reversed (Fig. 1(b)). All 
couplings goes to zero in the infrared and the interaction 
is irrelevant. The six Majorana fermions are massless and 
the model is critical with the central charge c = 6 x 1/2 = 
3. However, for generic values of the initial conditions 
(d(0) 7^ 0), the fixed point Hamiltonian does not have the 
SO (6) symmetry. Indeed, the velocities in both spin and 
orbital sector have different fixed point values m* ^ m^ so 
that the symmetry at the fixed point is SO(3)s(8) S0(3)t. 



It is only on the symmetric line (Gi = G2) that the fixed 
point symmetry is SO (6). 

- The C Phase. There one has /x(0) < and s can have 
both signs. As seen in Fig. 1(c), X{t) ^^ 00 when both 
t -^ ±cxD. The theory is not asymptotically free neither in 
the infrared nor in the ultraviolet. Though one certainly 
expects that a gap opens in the spectrum, perturbation 
theory is pathological [g4|. Indeed, from the field theo- 
retical point of view, the lack of asymptotic freedom in 
the ultraviolet implies that a well defined renormalized 
theory with three different couplings Gm, G2R and G^n 
does not exist. The low energy physics in such a situation 
is a highly non trivial and essentially non perturbative 
problem and one is left to make a sensible conjecture. 

As in the A phase discussed above, in the far infrared 
Gi(i) ~ G2{t) ^ — 2G3(i) as well as u* = Uj so one may 
wonder again whether the SO (6) symmetry is approxi- 
matively restored or not. We stress that the situation 
at hand here is very different to what happens in the 
A phase since the theory is not asymptotically free in 
the ultraviolet. This may have important consequences 
for the low energy physics. Indeed it is well known that 
the divergency of some couplings at high energy is rem- 
iniscent of the fact that some part of the interaction is 
irrelevant at low energy. Of course, perturbation theory 
alone cannot tell us which part of the interaction is irrel- 
evant and at present all what we have at hand to make 
a reasonable hypothesis is our perturbative results. The 
simplest scenario is that the SO (6) symmetry is approx- 
imatively restored provided the initial conditions are not 
too anisotropic (see however subsection C). As a support 
of this conjecture let us mention that this is what hap- 
pens in the C phase of the XXZ model where the Bethe 
ansatz solution tells us that the SU(2) symmetry is re- 
stored up to exponentially small corrections [^ . We 
are of course aware that our hypothesis is highly ques- 
tionable but it is the simplest one and could in principle 
be tested either numerically or experimentally. Indeed, 
the immediate consequence of the possible SO (6) restora- 
tion is that the effective low energy effective theory of the 
Hamiltonian (]^ in the C phase is approximatively that 
of the SO (6) GN model: 
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(40) 



where the effective coupling G is positive. The model 
( pO| ) is integrable. Its spectrum is known and consists of 
the fundamental fermion, with mass M, together with a 
kink of mass m = A//\/2 p^. At this point the question 
that naturally arises is how this enlarged SO (6) symme- 
try reflects in the spin and orbital correlation functions. 
The answer to this important question requires the com- 
putation of the exact dynamical correlation functions in 
the SO (6) restored massive phase. This could be accom- 
plished in principle by the form factors approach as in 



the frustrated two leg ladder considered in Ref. ( p5[). 
However this task is even more for the SO (6) case and 
goes beyond the scope of this paper. 



B. Phase diagram 

Let us now sum up our results and present the phase 
diagram associated with Hamiltonian (24). As even at 
the SU(4) symmetric point the coupling G3 is not zero, 
the best way to visualize the phase diagram is to fix the 
value of G3 and to look in the plane (Gi , G2) . As one can 
see from Fig. 2, there are two curves separating the three 
regions A, B and C which are given by the equation: 



/I — G1G2 + G1G3 -|- G2G3 



0. 



(41) 



In the region B all models are critical with approximate 
SO (6) symmetry The spectrum consists into six mass- 
less Majorana fermions with different velocities in both 
spin and orbital sectors. As one crosses the critical line E 
one enters in a fully massive phase (C phase) with an ap- 
proximate SO (6) symmetry. For large positive values of 
the couplings Gi and G2 one may finally enter the region 
A. Notice that there is no phase transition between the 
A phase and the C phase since both phases are massive 
but rather a smooth cross over. 




FIG. 2. Phase diagram for anisotropic couplings in the 
plane (Gi, G2) at a fixed value of G3. The massless phase 
B is separated from the massive C phase by the critical line 
E. The special point S on the symmetric line which is at the 
border between the C and the A phases is S0(6) symmetric 
with Gi = G2 = -2 G3. We plot also the two trajectories 
labelled 71 and 72 for which the scaling of the mass gap have 
two different qualitative behaviors. 



The most important feature of our phase diagram is 
that the two regions with G1G2 < which would have 
been massless in the absence of G3 do not survive an ar- 
bitrarily small value of G3. Therefore the c — 3/2 phases 
discussed in the previous section are not stable. Though 
the whole region is essentially massive there are still room 
for criticality since the B phase extends in both quadrans 
G1G2 < 0. From Eq. (41) the width of the critical re- 
gion is of the order IGaj. Therefore the main effect of 
G3 is to drive the system either to a fully massive phase 
(phase C) or to a fully massless phase (phase B), both 
with an approximate SO (6) symmetry. Of course, the 
phase diagram as depicted in the Fig. 2 is stricly valid in 
the small G limit. In particular higher order corrections 
may modify the shape of the critical curve E. We stress 
however, that its behavior in the vicinity of the SU(4) 
point at Gi = G2 = is given by our one loop result. In 
particular the fact that E crosses the symmetric line with 
a right angle will not change as one includes higher order 
in perturbation theory. However, as one goes to large 
deviations from the SU(4) point our effective theory will 
not apply since for large enough positive (respectively 
negative) Gi and negative (respectively positive) G2 the 
orbital (respectively spin) degrees of freedom order fer- 
romagnetically p,M. 



C. Effect of the anisotropy 

At this point one may wonder whether the anisotropy 
has no effect at all in the physical quantities. Apart from 
the velocity anisotropy in the B phase, one may still ex- 
pect some non trivial effect of the anisotropy in the C 
phase. Indeed, after all eventhough both c = 3/2 phases 
are unstable there should be some significant signature 
of the presence of the S0(3)i fixed point in the scahng 
of the physical quantities and in the finite size scaling. 
The very reason for this is that the SO (6) symmetry is 
restored dynamically with help of a marginal operator. 

-Cross-over and finite size scaling. Let us look at the 
RG flow in more details. The regions of interest are those 
with G1G2 < in which, as discussed in the previous 
section, if not for G3 one of the coupling would have been 
irrelevant. Then either the spin or the orbital degrees of 
freedom would have remained massless. In the following 
we shall concentrate on the orbital degrees of freedom. 
All the results that will be given can be straightforwardly 
extended to the spin degrees of freedom. 

Consider first initial conditions slightly above the 
phase B in the lower right quadran with G1G2 < (see 
Fig. 2) . We are then in the phase C with Gi (0) > and 
G2(0) < 0. As seen in the Fig. 3(a), the flow can be 
splitted into two "time" slices. 





FIG. 3. Qualitative RG evolution of the coupling con- 
stants in regions where 01(0)6*2(0) < and where the 
cross-over "time" defined by G2(to2) ~ exists. 



At first jGal strongly decreases and remains almost 
constant. Then, both spin and orbital degrees of free- 
dom weakly interact since |G(i 2)| >> IG3I. In the mean 
time, the coupling G2 increases and changes its sign at 
a time io2 where it vanishes. At "times" t < io2 the or- 
bital degrees of freedom do not know yet they shall enter 
in a massive phase and the system is on the influence 
of the SO (3)1 fixed point. Finally, as t -^ Im all the 
couplings blow up and one enters in the strong coupling 
region. The physical interpretation of this phenomenon 
is clear; it is the spin degrees of freedom that drive the 
orbital degrees of freedom away from criticality and it 
takes roughly exp(io2) RG iterations for the orbital sec- 
tor to escape from the bassin of attraction of the S0(3)i 
fixed point. 

The other interesting region is the portion of the B 
phase delimited by the critical curve S and the G2-axis 
in upper left quadran of Fig. 2. There Gi(0) < and 
G2(0) > 0. As seen in Fig. 3(b), it is the spin degrees 
of freedom that drive the orbital ones to criticality. The 
coupling G2 changes its sign and vanishes at a time to2 
and goes to zero in the limit t — )■ 00. In both cases, the 
cross-over "time" to2 is given by the implicit equation: 



X{t02) = 



m 



d(0)G3(0)- 



(42) 



Of course a similar discussion holds for the coupling Gi 
and defines another cross-over time ioi which is defined 
by the implicit equation: 



X{toi) = 



d(0)G3(0)- 



(43) 



We plot in Fig. 4 the iso-to2 and iso-toicurves in the 
plane (Gi,G2). 



Gi(toi)=0 




Go 



Gi 




G,(to2)=0 



FIG. 4. Iso-to curves in the plane (Gi, G2); the dashed hne 
corresponds to an iso-f 02 associated with the orbital degrees of 
freedom whereas the dashed doted line is an iso-toi associated 
with the spin degrees of freedom. E denotes the critical line. 




FIG. 5. Isogap curves in the plane (Gi,G2) 
the critical line. 



E denotes 



We stress that this cross-over phenomenon may have 
important practical consequences in numerical simula- 
tions. Indeed, in a finite system of size L, the criti- 
cal region will seem to extend towards the iso-ioi with 
toi = InL in the region Gi < 0, 6*2 > and towards the 
iso-to2 with to2 = InL in the region Gi > 0, G2 < 0. The 
low lying spectrum will look like if either the spin or the 
orbital degrees of freedom would have been massless. It 
is important to notice that these two pseudo-critical lines 
bend upwards in contrast with the critical line S. 

- The mass gap. The anisotropy has also non trivial 
effect on the mass gap of the system. To see this we have 
computed the mass gap M with help of the RG equations 
(p8|). As well known, the gap M is defined as the scale 
where perturbation theory breaks down. More precisely 
it is given by the scale tM = ln(l/Mao) at which the 
couplings blow up: 



X{tM) = 00. 



(44) 



The above implicit equation can be solved but the result- 
ing analytical expression is too cumbersome to be quoted 
here. Details are given in the Appendix and we shall 
content ourselves by its asymptotic behaviors as one ap- 
proaches the critical surface E. 



We find two qualitatively different behaviors depend- 
ing on the way one approaches S. When one approaches 
E at the SU(4) point (Gi = G2 = 0) we find: 



M~ Aexp(-G(7i)/A2/3) 



(45) 



where A = \/Gf + G\ -^ and the constant G(7i) de- 
pends on the trajectory labelled 71 in the Fig. 2. On the 
other hand, as one approaches S at any other point the 
asymptotic behavior of the gap is different: 



Af-Aexp(-G(72)/A) 



(46) 



where now A is the Euclidean distance of the critical sur- 
face. As above G(72) depends on the curve labelled 72 
in Fig. 2. There are two other marginal cases as one ap- 
proaches E tangentially: the exponents of A in both ( |45| ) 
and (Hq) are doubled. We see that the gap is generically 
larger in the regions with small anisotropy. To visual- 
ize this phenomenon we present in the Fig. 5 the isogap 
curves in the plane (Gi, G2). We observe that at a given 
distance A from the critical line E the gap increases as 
one moves towards the symmetric region (Gi ^ G2) and 
is maximum on the symmetric line (Gi = G2). Similarly, 
to keep the value of the gap to some constant M one has 
to move away from E as one leaves the symmetric region. 
-Effect of anisotropy on the S0(6) symmetry restora- 
tion. We conclude this section by looking at the effect 
of the anisotropy on the SO (6) symmetry restoration in 
the C phase. As discussed above, in the C phase the cou- 
plings Gi and G2 tend to — 2G3 in the far infrared (i.e. 
as i — *■ tjii) so that the effective Hamiltonian has an ap- 
parent SO (6) symmetry. Looking in more details at the 
asymptotic behavior of the couplings Gi [t) and G2 {t) we 
find that: 

G(i,2)(0 - -2Gmx\t) ± ^Xit). (47) 



Therefore, as when t -^ tM, X{t) -^ oo, there remains 
a subdominant infrared singularity when d(0) ^ 0. We 
thus expect some corrections to the SO (6) behavior. To 
get an estimate of these corrections we make the rea- 
sonable hypothesis that, when the anisotropy is small 
enough, the effective theory will be given by an SO (6) 
GN model with a small symmetry-breaking perturbation 
proportional to: 
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+ X ({k3 -I- K4 -I- K^f - (ki -I- K2 -I- Ke)^] , (48) 

where X << G. Of course we do not know the value of 
the effective coupling A, but it is sufficient for our purpose 
that it is small enough which will be the case if d(0) is 
small. From Eq. (48) one can get an estimate of the 



mass splitting in both the spin and orbital sectors: 



Mt 



1 



A 



In (l/Afflo) 



(49) 



Therefore, we expect the SO (6) symmetry will hold ap- 
proximatively only for small enough gap and anisotropy. 
The above argument is of course far from being rigorous 
but it helps to get an idea of the effect of the anisotropy 
in the symmetry restoration process. This should serve 
as a warning that, once again, one should be careful! with 
the conclusions drawn from perturbation theory. 



the spin degrees of freedom align parallel to the field and 
the remaining degrees of freedom will decouple. More in- 
teresting is the situation when h is small. In this case, one 
may expand ( pl| ) around the SU(4)i fixed point, as in pre- 
vious section, and study its effects in perturbation. For 
small values of h, the effective low energy Hamiltonian 
associated with (pw can be easily derived using the ex- 
pression of the spin operator S^ in terms of the Majorana 
fermions as given by Eqs. (p^),(pO|) and (21). Dropping 
all oscillating terms, the only contribution comes from 
the uniform part of the spin density and the effect of the 
magnetic field is just to add to the Hamiltonian ( |24| ) the 
interaction Hz = ih{^}^(^j^ + Ci'Ci) and the low-energy 
Hamiltonian thus writes: 



H/i = -^ UsB. ■ dxisB. ~ ^sL ■ dx^sLJ + 9l («! "f K2 + Kg)^ 

iXtR ■ dxitR - itL ■ dx^tLj + 92 (k3 + K4 + Kzf 



lUt 



+ 2Gz (ki + K2 + Kg) (k3 + K4 + K5) 



(52) 



where we recall that 



gi=Gi + G3 
52 = G2 + G3. 



A. Renormalization group analysis 



(53) 



IV. EFFECT OF A MAGNETIC FIELD 

In contrast with the SU(2) Heisenberg chain there is 
no unique way to apply a magnetic field in the SU(4) 
model. Indeed, while in the SU(2) case there is, apart 
from the total spin, only one conserved charges S^ in the 
SU(4) model there are three of them. These are the three 
Cartan generators of the SU(4) Lie algebra. In this work 
we have chosen for the three commuting generators: 



H^ = S\ H^ 



T^ and H-^ 



2gZrpZ 



(50) 



Other choices are possible but (|5y) is the one that is 
physically relevant to our problem. We thus see that one 
is at liberty to apply a magnetic field in any "direction" 
in the Cartan basis (pO|). However, if one relies on the 
physical interpretation of both operators S and T as spin 
and orbital operators, a magnetic field only couples to 
H^ — S^\ The resulting Hamiltonian thus writes: 

i 

+k(5, •5,+i)(f,.r,+i)+/z^5f. (51) 

i 

The magnetic field breaks the SU(4) symmetry down 
to SU(2)(g) U(l). For a sufficiently large value oi h > ho 



The situation at hand is similar to that of the XXZ 
model in a field |g6[. There exists a magnetic length A/j 
below which the magnetic field has no substential effect 
on the physics. However as the scale A >> A/j, the field 
strongly modifies the low energy behavior. Indeed, in this 
regime, as we shall see below, the two Majorana fermions 
C^ and C^ completely decouple from the four others. The 
best way to see this is to bosonize the "1-2" sector of 
the Hamiltonian (|5^) with help of Eq. (Q . The resulting 
Hamiltonian becomes: 



'^h = ^(^ (dx^sf + K {d^esf 



We 



&.eR - eLd.a) 



TT ( 6fl ■ dx^tR — itL ■ dx^tL 

+ 2/2 (k3'«4 + K3K5 + K4K5) 

+ 2/3 Ke (k3 + K4 + K5) 



2il 



nao 



cos f Vivr^s + 2qx] (K3 + K4 + K5) 



cos I V 47r$s + 2qx ) kq 



(54) 



where 
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TT 

Kh 
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(55) 



In the above Hamiltonian, we have introduced new veloc- 
ities and couphng constants with bare values: ug = Ws, 
vt = ut, f2 ^ 92, fs = fi ^ Gs and /s = gi correspond- 
ing to the operators that are stable under renormaliza- 
tion. As readily seen, the cosine terms involving $s will 
start to oscillate with wavevector q ^ h. Therefore, at 
sufficiently large magnetic field the two last terms in JS^ ) 
can be dropped and the bosonic field $s completely de- 
couples. When a gap is present, i.e. when on sits in the 
C phase a.t h — 0, there exists a critical field he ^ M 
below which the field has little effects on the low energy 
physics. As h increases and becomes greater than he, the 
spin sector undergoes a commensurate-incommensurate 
transition of the JNPT type |M to an incommensurate 
massless phase with central charge Cg = 1- In contrast, 
in the absence of a gap, i.e. when one sits in the B phase 
at h = 0, the spin sector always decouples at sufficiently 
low energies. Of course, this decoupling procedure has 
to be understood in the framework of the renormaliza- 
tion group. To this end we have computed the one loop 
recursion relation associated with the Hamiltonian (54). 
Neglecting the velocity renormalization we obtain: 

X = -(3/| + /|)Jo(2g) 

f2^f^ + f! + 2fiM2q) 

/3 = 2/2/3 + 2/4/5 Jo(2g) 

/4 = (1 - K)U + (2/2/4 + /3/5)(l + M2q))/2 

/s = (1 - K)h + 3/3/4(1 + Jo(2g))/2 (56) 

where the couphngs have been rescaled as / ^ J/t^v*, 
V* being a velocity scale, and J(o.i)(a;) are the Bessel 
functions of the ffi'st and second kinds. Apart from the 
inherent anisotropy in our problem, the above equations 
resemble to that of the XXZ model in a field |Q . We see 
on Eqs. ( |5^ ) that when q{t) ^ n the two Bessel functions 
start to oscillate and the renormalization comming from 
both the /4 and /5-terms is stopped. This defines the 
magnetic length Xh as q{th) = tt (th — InA/i) where the 
spin degrees of freedom decouple from the rest of the in- 
teraction (as well known, there is no unique way to define 
the magnetic length but we have checked that other rea- 
sonable choices do not affect qualitatively the physics). 
At this point, it is worth stressing that in the massive 
phase, th must be much smaller than Im, the scale at 
which perturbation theory breaks down. This means, of 
course, that h > M as discussed above. 

At times t > th, the couplings f^ and f^ do not partic- 
ipate in the RG equations of the couplings K, q, /2 and 
/s and Eqs.(p6h reduce to: 



and 





K = 




q = q 


/2 


-fl + fi 


/3 


= 2/2/3. 



(57) 



(58) 



This means that in the regime t > th the Hamiltonian 
(M) decouples: 



Ti-h ~Hs +'Ha 



(59) 



where 



w. = ^('-| id,^,f + K {d^e.f] (60) 



and 



-Ha 



ivt 



Kt-R • dx^tR — i,tL ■ dx^tL 

+ 2/2 (K3K4 + K3K5 + K4K5) 

-I- 2/3 Kg (^3 + K4 + K5) 



(61) 



where K,v^, wg, Vt, /2 and f^ are the effective couplings 
at the magnetic length A/j. 

Let us first concentrate on the spin sector as given by 
(|60|). The Hamiltonian ( pOJ ) is that of a Luttinger Liq- 
uid with stiffness K. The spin sector is thus massless 
and contribute to a central charge Cg — I. In addition, 
the correlation functions involving the field ^^ will be 
incommensurate with /i-dependent wavevector q. From 
Eqs. (21) one sees that the incommensurability will re- 
flect in both the spin and orbital correlation functions. 

We focus now on the remaining part of the interaction. 
The Hamiltonian 'H± describes the interaction between 
the orbital sector with the remaining spin-orbital Z2 de- 
gree of freedom (^^) of the spin sector. The low energy 
physics in this sector is non trivial and at issue is the 
behavior of the RG flow associated with Eq. (|5g) where 
the initial conditions have to be taken at the magnetic 
length with f2{th) and faith) at th = In(Ah). 

These equations are trivially solved. Indeed, upon in- 
troducing the new variables: 



/± = /2 ± fs 



(62) 



Eqs.(p8|) decouple: 



/± = /±- (63) 



As in the previous section, we distinguish between three 
phases A, B and C depending on the initial conditions of 
the flow: 

- The A Phase. This is when f+(th) > and f-(th) > 
0. Both couplings are relevant and a gap opens in the 
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spectrum. Moreover, since the theory is asymptotically 
free in the ultraviolet there are two length scales in the 
problem: m± -- exp-(7r//±(i,i)). 

- The B Phase. There f+{th) < and /_(i/i) < 0. The 
couplings are irrelevant and the four Majorana fermions 
become massless leaving the theory with the central 
charge c± = 2. The fixed point has only an approx- 
imate SO (4) symmetry since there remains a velocity 
anisotropy. Indeed, as in the zero field case, in general 
v^ ^ Vq . The generic symmetry of the fixed point is thus 
rather SO (3) (»Z2. 

- The C Phase. Finally is the C phase where /_(i/i) > 
and f+{th) < 0. Then f-{t) is relevant and f+{t) will go 
to zero in the infrared. Therefore, as in the previous sec- 
tion, one may conjecture that the SO (4) symmetry is ap- 
proximativelly restored. In the far infrared the effective 
Hamiltonian is that of the SO (4) Gross-Neveu model: 



Wi 



lU 



Y.^CRd.^R-CLd.CL) + f-{th) 



a— 3 



E 



txiitvj 



(64) 



which is integrable. Its spectrum consists only on kinks 
and antikinks with mass m ~ exp— (7r//_(</i)) plj . 

Notice, and this will be important for the discussion 
that will follow, that in both the phases B and C the 
effective theories are given (up to a velocity anisotropy) 
by (pj) with the difference that f-{th) is negative in the 
phase B (so that the interaction is irrelevant) while it is 
positive in the phase C. 



B. Phase diagram 

It is clear from the discussion given above that the val- 
ues of the effective couplings at the magnetic length are 
crucial. Of course, the existence of the commensurate- 
incommensurate transition in the spin sector depends 
only on the mass gap Af of the zero field Hamiltonian. 
What is more interesting, is what happens in the remain- 
ing orbital and spin-orbital sectors described by the Ma- 
jorana fermions 'C'^,C^,C^,C^- As we shall see now the 
anisotropy will play its tricks. Indeed, what is into ques- 
tion is the sign of the coupling f-{th) at the magnetic 
length. Returning to the original variables, one finds 
that /_ = G2. We saw in the preceeding section that 
in zero magnetic field the time evolution of G2 was very 
sensitive to the presence of the S0(3)i fixed point in the 
orbital sector and could change its sign at a time to2 de- 
pending on the the initial conditions (see Fig. 3(a) and 
Fig. 3(b)). The presence of a field h does not affect qual- 
itativelly this feature and provides for a renormalization 
of io2 which becomes /i-dependent: to2 — > tQ2{h). We 
consider now two cases. 

- First is when one sits in the B phase between the Gi 
axis and the critical surface S, with Gi < (see Fig. 
2). At zero magnetic field the system is critical with 



the central charge c = 3. There /_ is positive and de- 
creases as t increases vanishes at some RG time to2{h) 
and then changes sign. Now if t^ < to2{h), f^{th) will 
be still positive and a gap will open in the orbital and 
spin-orbital sector according to ([6^). On the other hand, 
if th > to2{h) then f-{th) < and there is no gap. This 
means that there exists a critical value of the field hco 
above which a gap opens. The portion of critical surface 
S in the region Gi < is thus unstable. The physical in- 
terpretation of this result is clear. In zero field, it was the 
spin degrees of freedom that drived the orbital degrees of 
freedom to criticality. When the field is large enough, its 
effect is to decouple a part of the spin degrees of freedom 
before the remaining fiuctuations had a chance to enter 
the bassin of attraction of the S0(4)i fixed point. Thus, 
the effect of the magnetic field in this region is to reduce 
the extension of the phase B. 

- The other interesting region, is when one sits in the 
C phase just above the critical surface S in the lower 
right quadran of Fig. 2, i.e. when Gi > and G2 < 0. 
There G2 is negative but is driven to positive values by 
the spin degrees of freedom. It changes sign at a time 
to2(/i) where it vanishes. Now if th < to2{h), f-{th) 
will be still negative while if th > to2{h), f-{th) will 
be positive. Therefore there exists a critical field /icO 
below which the orbital and spin-orbital sectors will be 
still massive. Above h^o, the gap will close. Again, the 
physical reason why the gap vanishes above h^o is that, 
the spin degrees of freedom did not have enough time to 
drive the orbital sector to strong coupling. We therefore 
conclude that the B phase has a tendency to extend in 
the region Gi > 0, G2 < when a field is present. 

To summerize, we expect two kinds of transition as one 
varies the magnetic field. When the theory is massive at 
/i = 0, as one increases h there will be a first transition 
in the spin sector to an incommensurate phase with the 
central charge c^ = 1. The remaining degrees of freedom 
will be still massive but are described by the SO (4) GN 
model. The coherent fermionic excitations of the SO (6) 
GN model disapear from the spectrum and the only mas- 
sive excitations that remain arc the kinks of the SO (4) 
GN model. Consequently all excitations will be incoher- 
ent. What happens as one increases h further strongly 
depends on the anisotropy. If G2 > the magnetic field 
just renormalizes the mass of the S0(4) kinks, the spec- 
trum is still massive. The total central charge of the 
model is thus c = Cs + c±_ ~ 1. However, when G2 < 0, 
there is a second phase transition at a field h^o of the 
Kosterlitz-Thouless (KT) type Eq] to a commensurate 
massless phase with cj^ = 2. Both spin and orbital de- 
grees of freedom are massless and the total central charge 
is c = Cs -I- c^ = 3. 

Notice that there will be three different velocities: Ug 7^ 
v*> ^ v^ so that the symmetry at the fixed point is not 
S0(6) but rather S0(3)(^ U(1)(^Z2. 

When there is no gap at zero field the spin sector is al- 
ways critical with incommensurate correlation functions. 
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What happens for the spin-orbital and orbital degrees of 
freedom depends again strongly on the anisotropy. When 
G2 < they remain massless and the total central charge 
is thus c = 3. However, if G2 > 0, there will be a KT 
type phase transition at a critical field hco to a massive 
phase with approximate SO (4) symmetry. 

We stress that the mechanism that leads to the KT 
type phase transition at the magnetic field hco is highly 
non trivial since the magnetic field does not couple di- 
rectly to both the orbital and the spin-orbital degrees of 
freedom. 



V. CONCLUSIONS 

In the present work we have studied the effect of sym- 
metry breaking perturbations in the one dimensional 
SU(4) spin-orbital model. Using the low energy effective 
field theory developped in Ref. ||lj| , we have investigated 
the phase diagram of the SU(2)(g)SU(2) model where the 
exchange in both the spin (Ji) and the orbital (J2) sec- 
tors are different. We found that the different phases of 
the symmetric Ji — J2 line extend to the case Ji ^ J2. 
In particular the massless phase, governed by the SO (6)1 
fixed point, extends to a finite region in the plane (Ji, J2) 
around the SU(4) point (Ji = -ft:/4, J2 = K/i). Simi- 
larly, in the vicinity of the critical surface, the massive 
phase has also an approximate S0(6) symmetry provided 
the anisotropy is not too large. In this phase, as in the 
isotropic case, the system spontaneously breaks transla- 
tional invariance and dimerizes with alternate spin and 
orbital singlets [|lj] . Both spin and orbital excitation are 
coherent at wave vector k = 7r/2. All these results remain 
valid in the vicinity of the SU(4) point. The question 
that naturally arises is what happens when K decreases. 
Indeed, in the limit K <^ "^(1,2) one enters in the weak 
coupling limit where magnon excitations are incoherent 
at wavevector k = tt [Q. In the simplest scenario, as dis- 
cussed in Ref. Q , one expects that the coherent peak at 
k = -K 12 in the dynamical susceptibility will disapear at 
a critical value oi K = Ku- Such a special point where 
an oscillating component of the correlation function dis- 
apear is a disorder point |29| ] and therefore, we do not 
expect a phase transition at K^ but rather a smooth 
cross over. 

Though these results could have been anticipated on 
the basis of the previous study of the symmetric case 
]14| , since the interactions are marginal, the anisotropy 
reveals itself in a non trivial scaling of the physical quan- 
tities. Indeed, we have shown that the anisotropy plays 
its tricks in two particular regions of the phase diagram 
with GiG'2 < 0, where G(\^2) = ^(1,2) — ^/4 measures 
the departure from the SU(4) point. In these regions, 
both spin and orbital degrees of freedom compete. For 
instance, when Gi > and G2 < 0, the spin sector tends 
to open a gap while the orbital one wants to remain mass- 
less. Since both sectors interact marginally, at issue is a 



delicate balance between the strength of the interactions: 
it is one kind of degrees of freedom that drives the other 
to its favorite behavior. In particular, this is the reason 
why the massless phase extends in the region with either 
Gi > or G2 > 0. In any cases, the whole system ul- 
timately becomes either fully massless or fully massive. 
The crucial point is that since the interactions are only 
marginal, it may take a very long time, in the renormal- 
ization group sense, to reach the asymptotic low energy 
regime. This has important consequences. 

First of all is the non trivial behavior of the mass gap 
of the system. We found that the gap M is genericaly 
smaller in the regions with large anisotropy, i.e. in the 
two quadrans G1G2 < above the critical curve S. This 
is due precisely to the strong tendency to massless be- 
havior in these regions. As a consequence the gap AI 
has two qualitatively different scaling behaviors as one 
approaches E either from the symmetric region or the 
asymmetric one (the trajectories labelled 71 and 72 in 
Fig. 2): M - exp(-Gi/A2/3) and M ~ exp(-G2/A). 

Second is the finite size scaling. Since the gap opens 
exponentially it is very difficult to localize accurately the 
critical line S in a finite system. In the current model the 
situation is even more ackward in the regions G1G2 < 0. 
Indeed, in a finite system of size L, the critical region will 
seem to extend and the pseudo-critical lines will be given 
by the two iso-ipi and iso-to2 curves, with ^0(1, 2) = InL, 
that have the opposite curvature than the true critical 
line S (see Fig. 4.). In this pseudo-critical region, ei- 
ther the spin or the orbital degrees of freedom will look 
massless. The phase diagram in zero magnetic field as 
obtained by very recent DMRG calculations ||l^,Q is in 
qualitative agreement with our RG analysis. However, 
the critical line obtained in these numerical computa- 
tions has the opposite curvature that the one loop result 
S. Our interpretation of this fact is that what has been 
observed are the two iso-io(i.2) curves. This reflects once 
again the non trivial finite size scaling induced by the 
anisotropy. 

Finally, is the effect of a magnetic field. The magnetic 
field affects the spin degrees of freedom in the usual way. 
In the massless phase it leads to incommensuration in the 
spin sector while when a gap is present, a commensurate- 
incommensurate transition can occur at a critical field. 
However, what happens to the remaining degrees of free- 
dom strongly depends on the anisotropy. In the region, 
where both degrees of freedom do not compete, i.e. when 
G1G2 > 0, the remaining orbital and spin-orbital sector 
remains either massless or massive with an approximate 
SO (4) symmetry. On the other hand, the most strick- 
ing effect occurs when the spin and orbital fluctuations 
compete i.e. for G1G2 < 0. In this region the held rein- 
forces the effect of the orbital degrees of freedom and can 
induce a second phase transition, of the KT type, for a 
sufficiently large field, from massive to massless approxi- 
mate S0(4) behavior. The origin of this non trivial effect 
of the magnetic field stems from the interplay of the pres- 
ence of orbital degeneracy and anisotropy. We hope that 
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this transition will be observed in further experiments on 
quasi one dimensional spin gapped materials with orbital 
degeneracy. 

Note added. After this work was completed, we became 
aware of a work by Itoi et al pO] who also predict the 
extension of the massless phase in the anisotropic region. 



VI. APPENDIX 

In this Appendix, we shall compute the mass gap M 
in the phase C, and obtain the asymptotics of M in 
the vicinity of the critical surface S. As well known, 
within perturbation theory the gap defines the scale tu — 
ln(l/Afao) where all the couplings blow up. Clearly, Im 
is given by the equation: 



X{tM) 



(65) 



Integrating ( |37| ) and recalling the dynamic of X in the 
phase C we find: 



^M — 



X* 



dX 



X^P(X) 



(66) 



where e — sign(s); P{X) is given by Eq.(|3q) and has only 
two reals roots and X* < 1. In the following s,fj. and 
G3 have to be understood as initial conditions. 
Performing the integrals we obtain: 



I'M 



[Eia{0),k)~eE{a{l),k)] 



lJ^\[j-^) [^("(0), fc) - eFia{l),k)] 



G3 



1 



P 



2|G3|(p- 



1 



(67) 



where F and E the elliptic functions of the first and the 
second kind respectively, with parameters: 



a{u) — 2 arctan 



Ip + u* - Gl/2^iu*^ 



2p 



(68) 



u* and p^ being given by: 

u* = 1/X* 



P 



2G| 

fJLU* 



(69) 



Developping Im arround a point (Gic, G2C, Gsc) be- 
longing to the critical surface E between the phases B 
and C (see Fig. 2), we obtain the asymptotics 

M-Aexp(-G(7i)/A2/3) ifG^i^G^^O 

M-Aexp-G(72)/A if G^=^G§ (70) 

where A is the Euclidian distance from S. The two con- 
stants G(7i) and G(72) are given by: 



G(7i)= 0.6845 G~//^ (cos 6*) -2/3 



G2c)') cos2 ( 



C(72) = ((l + 2G3cV(Gic 
where 9 is the angle to the normal of S 



-1/2 



(71) 
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